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"Take the world from another point of view." - R. Feynman 



Abstract. The oblique basis method is reviewed from engineering point of view related to 
vibration and control theory. Examples are used to demonstrate and relate the oblique basis in 
nuclear physics to the equivalent mathematical problems in vibration theory. The mathematical 
techniques, such as principal coordinates and root locus, used by vibration and control theory 
engineers are shown to be relevant to the Richardson - Gaudin pairing-like problems in nuclear 
physics. 



1. Introduction 

It is well known and commonly experienced observation that "... mathematical concepts turn 
up in entirely unexpected connections ..." The surprise of finding a mathematical structure 
applicable to seemingly unrelated fields of study is somehow proportional to our ignorance of 
one or the other subject, even though we may have studied the underlining mathematical field 
previously. The following contribution attests to this situation by drawing the analogies and 
pointing examples from my research topic in nuclear physics and a more recent engineering course 
in vibration and control theory which I had the pleasure to teach. This reinforces my view that 
one should keep an open mind and must often try to "Take the world from another point of 
view" [2]. In this respect, this article is devoted to the discussion of two major applications 
of a linear algebra: matrix diagonalization and location of zeros of complex algebraic systems. 
The first application is in the subject of complex many-body systems like nuclei ^ while the 
second is their application in the filed of vibration and control theory [4] . Although presenting 
the classical mechanics problem followed by its quantum mechanics framework within which the 
mathematical tools reside may be the better pedagogical point of view, we will use the opposite 
order of the exposition, due to the fact that the audience to which the material was presented 
is more familiar with my nuclear physics research framework. 

In what follows, the article justifies the oblique basis approach to complex problems and 
introduces the two important collective models commonly used in nuclear physics. Then the 
oblique basis method is stated as it is pertained to nuclear physics. The engineering problem 
is stated next and the relevant mathematical techniques, such as Evans root-locus method [5] 
for studying the stability of a system, are outlined. Finally, we go back to the nuclear physics 
subject and show the relevance of the Evans root-locus method for the study of exactly solvable 
nuclear pairing interaction. 



Before continuing the article with its main subject, I would like to acknowledge that, even 
though the oblique basis idea has been developed mostly as a collaborative effort O El [7], in 
my view it has been stimulated and supported strongly by my Ph.D. adviser Professor Jerry 
Draayer, in whose honor this conference has been held. 

2. Important Exactly Solvable 2-body Nuclear Interactions 

Exactly solvable problems such as harmonic oscillator and other group-symmetry based models 
in physics have been essential in the qualitative and quantitative understanding of physics 
phenomena. The need of better, more accurate, and easily replicable quantitative results has 
naturally introduced the use of powerful computers to the physics research. Never the less, 
exactly solvable problems are still essential to our understanding of nuclear structure [8] and 
reactions [7], especially in the limit of weak perturbation to the exact limit [6j. In particular, 
the nuclear shell-model in the limit of independent single particle motion is exactly solvable; 
however, when a two-body perturbation interaction is added then it becomes the cornerstone 
of the configuration mixing approach to the corresponding many-body problem [8]. When 
the interaction between the nucleons becomes significant the perturbation method is no-longer 
applicable and cannot reproduce the correct behavior of the system. Examples of such cases 
are well-deformed nuclei where the rotational structure of the spectra becomes more prominent 
or near-closed-shell nuclei where pairing phenomenon take over the low-energy spectra of the 
nuclei. 

By using Fock space techniques within the independent particle model with traditional single- 
particle level labels, one can consider two important collective interactions in nuclei. These 
interactions represent two well known exactly solvable systems: 

• the Elliott's SU(3) model with degenerate single-particle harmonic oscillator energies 
(ej = Ej = e) and quadrupole-quadrupole interaction Q ■ Q [8j: 

H = eN -^Q-Q, 

• the nuclear pairing (P'^P) interaction with non-degenerate single-particle energies (ej): 

H = EiNi + GP^P. 

3. T=l Proton-Neutron Pairing - SO(5) RG-model 

While the exact solvability of the Elliott's SU(3) model is more or less well understood from 
group-theoretical point of view, the exact solvability of the particle number conserving pairing 
interaction has much less known implementations, even though the solution has been published 
since 1962 by Richardson [9]. The exploration and implementation of the Richardson- Gaudin 
particle number conserving pairing interaction may have been hindered by its more complicated 
algebraic structure as seen in the work of Gaudin [lOj and Ushveridize [11]. However, the 
advancement of the computer power and recent understanding of the relevant mathematical 
and computational structure has brought a renewed interest in the topic [12^ [T3] along with the 
understanding of the integrability of the pairing interaction [14[ [T^ . The pairing has even been 
extended to an A-body generalized pairing interaction [16J. Finally, the full T=l proton-neutron 
pairing has been formulated within the frame work of the Richardson-Gaudin algebras |17j . 

Here we briefly summarize the main results on the Richardson-Gaudin pairing with the 
purpose of pointing out some of the possible practical implementation issues. We recall the 10 
generators of the S0(5) algebra in a representation suited for nuclear physics problems [TT]. 
The three T = 1 pair-creation operators are defined as: bLu = nlfij, Lq- = {h\p\ + p\n\) / ^/2, 
and 6^]^ j = p\p\, where n and p refer to neutrons and protons, respectively, and i labels a 



single-particle basis (with i its time-reversed state) that may be associated with the spherical 
shell-model basis i = jm or with an axially-symmetric deformed basis i = am. The three pair- 
annihilation operators are 60, and The three components of the isospin operator 
[r+,i = {plhi+plhj)/y/2, Tb,i = {■PiPi+ptpt)/2 - {h\hi + h\hi)/2, and f_^i = {h\pi + h\pi) / ^/2] 
close the SUr(2) subalgebra of S0(5). These 9 operators together with the number operator 
= V^iPi + p\vi + + n\hi define the S0(5) algebra. The pn-pairing Hamiltonian with a 
pairing strength g and isospin breaking A is given by the expression [17j : 

H = Y.^,{N,+/\N,,) + ^-f-f + gY^bl^J,,,^,, (1) 



with the eigenvalues of ([T]) given by the expression: 



a=l 13=1 j 



^{2 + A) -At, 



+ |ro(ro - 1), (2) 



where Vi is the seniority of each i level, i.e. the number of fermions not paired in time-reversed 
states with isospin T = 1, tj is the isospin of the unpaired fermions, t = tj, M is the number of 
T = 1 time-reversed pairs, and Tq is the z-axes component of the total isospin, i.e. the eigenvalue 
of the operator Tq = J2i To,i- The total number of nucleons is = Np+Nn = 2M+J2i '^j whereas 
their difference is Np — = 2Tq. The quantum numbers M, Tq, Vi, and tj are conserved; T is 
also conserved if A = 0. 

The variables Cq and lj^ are solutions of the corresponding Richardson-Gaudin equations [17J: 

1 ^ lj^ v. + 2U-2 ^ ^ 2 1 

9 2 .^^ Zi - ea Q,/(^o)=i ~ 13=1 ~ ^" 

A M + To+t r, Mr, L r,, 

_ = Y Z Y Z y (4) 

where variables Zj are related to the single particle energies £j (zj = 2£j ) and the indexes are 
specialized to a spherical basis i = jm. 

Thus, the problem of solving the pairing hamiltonian is reduced to solving the above non- 
linear equations for Cq and W/j. In practical terms, without computers in the old days, this 
would mean hard work by hand. Nowadays, there are sufficient computer power to solve these 
equations by starting at a small values of coupling constants g and A and gradually increasing 
these values to the desired magnitudes. However, there is an apparent singularity problem when 
some of the pair energies become close to the a single particle energy level, that is, when some of 
the denominators get close to zero. This problem can be resolved by trial and error search on the 
values pass the singularity point - this is similar to the problem of finding initial configurations 
for the values of and cj/j to start the process with. Alternatively, one can use small complex 
component of the single particle energies or the Evans root-locus method to be touched upon in 
the last section. 



4. Oblique Basis in Nuclear Physics 

In second quantized form the effective Hamiltonian for a nuclear system is: 

H = Y^iO-to'i + ^ Y (5) 



Here, £i are single-particle energies derived from excitation spectra of one valence particle system 
and VkiAj are two-body matrix elements viewed as a residual interaction beyond the single- 
particle mean-field confining potential. 

When the residual interaction can be treated as perturbation then one can apply perturbation 
theory using the basis of the exactly solvable limit. If there are two or more exactly solvable 
limits then there are naturally two or more possible perturbation treatments that are based 
on their corresponding basis sets. However, if the strengths of the interactions are such that 
perturbation theory is not applicable, that is, there is a strong mixing, then a combination of 
relevant basis sets can be used to approach the problem. Thus, the oblique basis consists of two 
or more basis sets. In this sense the oblique basis idea has been applied to the case when the 
Elliott's SU(3) is broken due to the non-degenerate single-particle energies [3j and also illustrated 
on simple toy model of a particle in a box [6j. 

In our simple example [3] the two basis sets were: spherical shell-model states (ssm-states) 
expressed in spherical single-particle coordinates (n/j) that are eigenvectors of the single particle 
non-interaction hamiltonian, while the second set has a good SU(3) structure (su3-states) that 
track nuclear deformation. The SU(3) basis set is given in cylindrical single-particle coordinates. 
Such oblique calculations result in a nonorthogonal oblique basis and require an evaluation of 
the matrix elements of physical operators plus a knowledge of the scalar product (cq, - Ei) related 
to the overlap matrix. Schematically, these basis vectors are labeled by a = 1,..., dim(ssm-basis) 
and i = 1,..., dim(su3-basis) and the corresponding matrixes can be represented in the following 
way: 



, . , f Sa '■ ssm — basis \ 

basis vectors : ^ o i ■ K (o 

\ Ei : su3 — basis / 

overlap matrix : = ^ ^ ^ , = ■ Ei, (7) 



Hamiltonian matrix : H=( ^^««'™x^«™ ^^ssmxsu^ ^ ^^ap -aj _ 

In general the oblique basis results in a non-orthogonal basis that may even become 
linearly dependent and over-complete. Thus, the matrix equation includes an overlap matrix 
(0jj = and has the form: 

^ {H.jVj - \Q,jVj) = 0. (9) 
j 

In analogy to the SU(3) symmetry breaking in the pf-shell, which provides symmetry-adapted 
guidance to the relevant collective basis, one can expect that in a multi-shell space one would 
naturally use the symplectic symmetry group [THl US] . 

Furthermore, the oblique basis method can be extended and applied to the study of particle 
transfer reactions in strongly deformed nuclei where harmonic oscillator bound states techniques 
can be employed to compute the single particle binding energies in mean-field Woods-Saxon 
potential and then the proper tail of the wave function can be constructed using Sturmian basis 
that results in an equivalent oblique basis [7j. 



5. The Engineering Problems in Vibrations and Control Theory 

A system in mechanical engineering is often represented by a set of masses, springs, and dumpers 
[1]. The state of the system corresponds to a vector x that describes the coordinates of the 
masses. For our simple example we can neglect the dumper effects which usually corresponds 



to velocity dependent force [c]x. In this case the Newton's second law of motion is written as a 
matrix equation using mass [m] and spring [k] matrices: 

[k]x + [m]x = 

5. 1 . Normal Modes for a Simple System 

By using the method of separation of variables and the reasonable assumption that the system 
will exhibit some periodic motion near its equilibrium position, one sets x = XT{t) and 
T = —LtP'T. The resulting equation for the frequency of oscillations uj becomes effectively a 
generalized eigenvalue equation when the mass matrix is not diagonal. Simple example of such 
non-diagonal [m] is a wagon with a compound pendulum inside [4j. 
Solutions and vectors Xi to the generalized eigenvalue equation: 

\k\Xi-uj1\m\Xi = ^, (10) 

form a [A;] or [m] "orthonormal basis" when ujj ^ uji {Xi [m]Xj = and Xi [k]Xj = 0). Thus, 
in terms of the normal mood matrix X = {Xi, X2, . . . , X„} the super-matrices: M = X'^[m]X 
and K = X'^[k]X become diagonal. Since the mass matrix is positive definite due to its 
relation to the kinetic energy of the system, then one can consider [m] normalized normal 

modes Xi [m\Xj = 5ij which results in a diagonal super-matrices Mij = SijLof. In terms of the 
normal modes the system of n-coupled masses becomes a system of n-uncoupled masses which 
is equivalent to the study of n one-dimensional oscillators. Identifying the normal modes of a 
mechanical system is very useful once one introduces back general damping and other possible 
forces to the problem. Such additional forces can be treated in a perturbative like regime using 
the corresponding normal modes of the system. 

5.2. Evans Root-Locus Method 

How to influence the behavior of particular system to achieve the desired response is the subject 
of the Control Theory. In particular, by going back to our simple system but keeping only one 
degree of freedom we have a simple second-order differential equation. From the general theory 
of these equations, we know that if the system is displaced from its equilibrium there would be 
some oscillations and some damping that will determine the transition behavior. The response 
of the system to a Dirac delta input is the transfer function. Since any function can be viewed 
as a collection of Dirac-delta impulses the output to a general input is essentially related to the 
structure of the transfer function which is in essence the ratio of the output to the input signal 
but in the Laplace transform variable s. Usually systems have a proportional controller K and 
a response of the controlled plant G(s) that is a simple polynomial function of s which result in 
an overall transfer function T{s) = KG{s). 

A system is stable if all the poles of the transfer function l'{s) are in the left half complex 
plane. A pole that has a real negative value corresponds to a pure damping without oscillations 
while a pole with non zero imaginary part corresponds to an oscillatory behavior. If the system 
has a pole with a positive real value then the system become unstable and the amplitude of a 
small input signal can grow to infinity. To prevent such behavior many systems include sensors 
H{s) and negative feedback loops that reduce such output grow. Including a negative feedback 
modifies the transfer function of the system to become 7~(s) = KG{s)/{l + KG{s)H[s)). 

Finding how the poles of the transfer function move in the complex plane is essential for 
controlling the response of a system. Often the combination of transfer functions of various 
control components results in a transfer function that is a ratios of polynomials. Were the 
controller has a simple proportional response K to the input signal. The poles of the transfer 



function T{s) = KG{s)/{l + KG{s)H{s)) are clearly related to the zeros of I + KG{s)H{s) = 0. 
Understanding the zeros as K changes is thus essential. If one sets G{s)H(s) = L{s) = then 
Evans root-locus method gives rules for the location of the roots of the equation: 

Hs) = 4 

as K grows from 1 to infinity. Since L{s) is considered to be a proper fraction the polynomials 
b{s) = Yli^ii^ — Zi) and a{s) = YVi=i{s — Pi), that is, the degree n of a{s) is bigger or equal to 
the degree m of b{s) then L{s) has the form: 

This equation the main equation of the stability analysis in control theory based on Evans 
root-locus method [S]. 

The above mathematical expression and its corresponding mathematical problem is 
practically the same as our Richardson-Gaudin pairing equations discussed earlier ([3]). In 
particular, if we restrict the pairing to the problem of identical particles. Then we get the 
Richardson original equations [9j that only contains the pair energies 6q- and the single particle 
energies Zi = 2ej as some of the poles of L: 



1 

9 



L 

1=1 
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M 

E 



(12) 



6. Conclusion 

In the preceding text we have discussed the use of non-orthogonal basis sets as they emerge in the 
studies of complex systems like atomic nuclei and mechanical systems. In studying the nuclei, 
the motivation is to deal with strongly coupled competing modes where the system is far from 
any of the possible perturbation regimes. The well-known quadrupole-quadrupole and pairing 
interactions are examples of collective modes that can compete against each other as well as 
against the single particle spin-orbit interaction which is often considered as the main source of 
lifting the single particle energy degeneracy. In a sense, each of these three exactly solvable limits 
represents a " normal" mode of the system and can be used to address the strong coupling regime 
of the system. The similarity between the principle normal modes of a mechanical system and 
the oblique basis in nuclear physics is further enhanced by the analogy between the Richardson- 
Gaudin pairing equations and the Evans root-locus method from the control theory for studying 
the satiability of a mechanical system. The author hopes that the reader will find the analogies 
interesting and this work will stimulate future collaborative research on these topics and their 
implementations to nuclear physics and engineering problems. 
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